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Motivated by recent experiments in optics and atomic physics, we derive an averaged nonlinear 
partial differential equation describing the dynamics of the complex field in a nonlinear Schrodinger 
model in the presence of a periodic nonlinearity and a periodically- varying dissipation coefficient. 
The incorporation of dissipation is motivated by experimental considerations. We test the numerical 
behavior of the derived averaged equation by comparing it to the original nonautonomous model in 
a prototypical case scenario and observe good agreement between the two. 



I. INTRODUCTION 

In the past few years, there has been an intense theo- 
retical interest in the use of nonlinear Schrodinger (NLS) 
equations to describe both thepropagation of optical 
beams in waveguides and fibers [l|, HJ and the mean-field 
evolution of Bose-Einstein condensates (BECs) [H, 0|. 
Within this framework of dispersive equations that sup- 
port solitary nonlinear waves, one of the particular top- 
ics of recent interest has been the effects of spatially 
and/or temporally (i.e., in the evolution variable) depen- 
dent nonlinearities. This subject, often called "nonlin- 
earity management" [H — by analogy with the topic of 
"dispersion management" that has been developed, in 
the context of optics, in far greater depth [|| — was orig- 
inally proposed in the study of layered optical media pj . 
However, it has also garnered considerable attention in 
the study of Bose-Einstein condensation, where it was 
reformulated as Feshbach resonance management Q. 

Recent experimental work in optics has realized lay- 
ered media through a concatenation of glass slides and 
air gaps. This has allowed a more detailed examination 
of topics such as the breathing of localized pulses M and 
the modulational instability of extended ones [T(J EH- 
Moreover, in the context of Bose-Einstein condensates, 
the inter-atomic interactions (which are the source of the 
nonlinearity at the mean-field level) can be adjusted ex- 
perimentally in a very broad range by employing either 
magnetic [H, [l3| or optical Feshbach resonances [lij. 
This has led to a significant number of both theoret- 
ical and experimental studies, including the formation 
(in the laboratory) of bright matter-wave solitons and 
soliton trains for 7 Li [H, fl6[ and 85 Rb [l7| atoms. On 
the theoretical side, such a modulation of the interac- 
tion scattering length (and hence of the nonlinearity co- 
efficient) has been used, among other things, to stabi- 
lize attractive higher-dimensional BECs against collapse 
[Til ]. More recently, spatial variations of the nonlinear- 
ity have also been considered. In particular, it has been 
shown that such "collisionally inhomogeneous" conden- 
sates lead to a variety of interesting features, including 
adiabatic compression of matter-waves 0, H3, Bloch 
oscillations of matter- wave solitons [l9j], atomic soliton 
emission and atom lasers 1211 , enhancement of transmit- 



tivity of matter- waves through barrie rs [22l . |23| , dynami- 
cal trapping of matter- wave solitons [22j, stable conden- 
sates exhibiting both attractive and repulsive interatomic 
interactions [24| , the derealization transition of matter 
waves [25| . and more. Numerous different types of spa- 
tial variations of the nonlinearity have now been con- 
sidered, including linear 0, [53], random [2r| periodic 
[H H3, HI, and localized (step-like) [HI iHUfl] ones. 
There have also been a number of detailed mathematical 
studies [U HH, HH that address aspects such as the effect 
of a "nonlinear lattice potential" (i.e., a spatially periodic 
nonlinearity) on the stability/instability of solitary waves 
and the interplay between drift and diffraction/blow- up 
instabilities. 

In the case of fast variations of the nonlinearity coeffi- 
cient as a function of the evolution variable (time in BECs 
and the propagation direction in optics), one successful 
strategy that has been employed is to average the nonau- 
tonomous, nonlinearity-managed dynamics to obtain an 
(averaged) autonomous system [3 4| . [ 35| . The stationary 



states [36j and collapse properties [37( of the latter can be 
analyzed in one- and in higher-dimensions, respectively. 
However, the presence of dissipation is a particularly im- 
portant feature that arises when periodically varying the 
nonlinearity coefficients in both optics and BEC exper- 
iments and which was not incorporated in these earlier 
studies, to the best of our knowledge. In particular, in 
the optical case of nonlinearity management there is a 
periodic loss (of the order of a few percent of the inten- 
sity of the optical beam) every time the beam crosses an 
interface between the different media (such as glass and 
air) 0, [ToL |Tll |. In BECs, if the Feshbach resonance is 
crossed, numerous atoms are lost, which results again in 
dissipative dynamics 0, 0] (although it is important to 
note that it is not always necessary to cross the actual 
resonance to change the sign of the scattering length, as 
vanishings of the scattering length as a function of the 
external magnetic field can also occur away from the res- 
onance d, EH). 

Our goal in this Brief Report is to address the averag- 
ing of nonlinearly-managed dynamics in the presence of 
dissipation. The remainder of our presentation is orga- 
nized as follows. We first present the general setting in 
which both a periodic nonlinearity and a periodic dissi- 
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pation are applied. We then use averaging techniques to 
this nonautonomous setting and obtain an autonomous 
partial differential equation (PDE) describing the aver- 
aged dynamics. We subsequently test the resulting model 
against numerical experiments of the original dynamical 
equations, obtaining good agreement between the two. 
Finally, we summarize our findings and present some sug- 
gestions for possible future studies. 



II. ANALYTICAL RESULTS 

Motivated by the above physical settings, we consider 
in our analysis a time-dependent nonlinearity and include 
a time-dependent dissipation term. We thereby general- 
ize the averaging technique of [H| and obtain a general 
averaged PDE in arbitrary dimensions. As our derivation 
hinges on the periodicity of the fast timescalc, we require 
that the length of the period is at least an order of mag- 
nitude smaller than that of the slow scale over which we 
monitor the dynamics. 

The primary model used for the physical settings we 
consider is an NLS equation of the form 

iu t + \Au + la \u\ 2 u + iCou + 77(f) M 2 « + 7C(f)" = 0, 

(1) 

where x G M™, t E R+, r = t/e and the quantities 70, 
Co are parameters. The continuous functions 7(7"), C( T ) 
satisfy 

7 (r + 1) = 7 (r), /„ 7 (r)dT = , 
C(r + l) = C(r), / o 1 C(r)rfr = 0. 

and represent the time-dependent part of the nonlinear- 
ity and dissipation, respectively. The NLS equation de- 
scribes the envelope of the electric field of light in the 
context of optics, and it represents the mean-field wave- 
function of the BEC in the context of atomic physics. 
To better understand the behavior of solutions, we use a 
multiple-scale expansion to derive an averaged equation 
for (fT]) in arbitrary dimensions. 

Following the notation in Ref. [34J, let /_i denote the 
zero-mean antiderivativc of /. It is given by 

/_l(r) = f /(r')rfr' - / f /(r')dr'dr . (2) 
Jo Jo Jo 

Define the transformation 

it(x, t) = e~ c - l(T) e fo ' (r)|t,(x ^ )|2) v(x, t) , (3) 

where j{r) = (je-^- 1 ^) _ y Using (|2|3]1 , equation {J) 
can be expressed as 

iiH —j(r)\v\tV = -\Av - 7 e~ 2C - l(r) H 2 v - i( v 
(t) [2 (VM 2 • V«) + v am 2 ] 
+ ±j(t) 2 v(V\v\ 2 -V\v\ 2 ) , (4) 



where \v | 2 = ^(M 2 ), A|v| 2 stands for A(|z;| 2 ), and V|w| 2 
stands for V(|w| 2 ). 

We isolate \v\ 2 by considering the expression i;(|3])— 
[i.e., Eq. (j4]) times the complex conjugate of v minus v 
times the complex conjugate of Eq. Q] to transform ([TJ) 
into the standard form 

iv t = -±Av--t e- 2C - l{T) \v\ 2 v-iCov-2( 0] (T)\v\ 2 v 
-\ij{r) [A(\v\ 2 v) - 2\v\ 2 Av + v 2 ~Av] 
-ij(r) 2 [2\v\ 2 A\v\ 2 + W\v\ 2 -V\v\ 2 ]v. (5) 

Using the multiple-scale expansion v(x,t) = w(x,t) + 
evi(x,t,r) + 0(e 2 ), we obtain from ([5]): 

iw t = — iv\ T — iAw — 7oe _2< " _1 ( r ^|w| 2 u' 
-i( w - 2( j(T)\w\ 2 w 
-\ij{T) [A(|w| 2 w) - 2|w| 2 Aw + w 2 Aw] 
-ij(r) 2 [2| W | 2 AM 2 + VM 2 • VM 2 ] w ■ (6) 

Integrating ([5]) yields an expression that one can reintro- 
duce into the equation to solve for v± T . Consequently, 
the averaged equation for |T]) takes the form 

iwt — -|Ato — r yop\w\ 2 w — iCo w 

-^-w [y\w\ 2 ■ V\w\ 2 + 2|w| 2 A|w| 2 ] , (7) 

with a 2 — j(r) 2 d,T and p = e^ 2 ^- 1 ^^. Observe 
that the formal expansion v = w + ev\ yields an equation 
which no longer depends upon the fast time-scale r (that 
is, an autonomous PDE). This approach also enables us 
to obtain the governing dynamics for the leading-order 
correction to the averaged behavior: 

vi = -\ [A(|w| 2 w) - 2\w\ 2 Aw + w 2 Aw] j-i(r) 

+ \iw [V\w\ 2 ■ VM 2 + 2H 2 AM 2 ] (j(T) 2 - a 2 )_ 1 

+i~fo\w\ 2 w (e" 2C - l(T) )_ i + 2iCoJ-i(t)\w\ 2 w , (8) 
which can be compared with Eq. (2.13) in Ref. [34| . 

III. NUMERICAL CORROBORATION 

In order to test the validity of Eq. for the averaged 
dynamics, we implement a prototypical two-dimensional 
realization of the above setting with radial symmetry (fol- 
lowing the lines of [HI). In particular, we consider the 
case of the two-dimensional unstable (against collapse) 
NLS soliton — the so-called Townes soliton [38| ) — with 
focusing nonlinearity. It has been shown in the context 
of BECs that such a solution can be stabilized using a 
rapidly oscillating nonlinearity coefficient [l8j]. The re- 
sulting dynamics contain a fast time-scale periodicity as- 
sociated with the nonlinearity management, so that this 
setting provides an ideal testbed for examining the accu- 
racy of Eq. 0. 
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FIG. 1: (Color online) Averaged dynamics of Eq. (0 versus 
full dynamics of Eq. the maximum of the intensity of 

(£9j) is shown in gray (green in the online version). We use 
the parameters ao = 2ir, ai = Sir, u> — 30, b = bo x 10~ 4 , 
and bo = 2. Overlayed (in black) is the same diagnostic for 
the averaged equation. It is clear that the averaged dynamical 
behavior is captured accurately over 250 time units. The inset 
shows a magnification of the small time-scale oscillations that 
are responsible for the stabilization and the overlayed average 
curve over the first 12 time units. 



As in Refs. [18J, we consider the setting 



FIG. 2: (Color online) Three comparisons like that of Fig. [T] 
over the first 60 time units with bo = 1,2, and 3 for the top, 
middle, and bottom curves, respectively. All other parame- 
ters are the same as in Fig.Q] In particular, observe in the top 
comparison that the averaged equation (0 begins diverging 
from the long time-scale dynamics (in the case of the smaller- 
magnitude dissipation term). The inset shows the absolute 
error between the solution of the full equation, Wfuii, and that 
of the averaged, u avc , Error= max(|u avc | 2 ) — max(|uf u ii| 2 )|(t). 
One can see more clearly here that after t = 30 the error in 
the case with the smallest dissipation coefficient grows con- 
siderably. 



id t + l -A +1 (t)\u\ 2 +iai) 



u = 0, 



(9) 



where the time-dependent nonlinearity 7(i) and the time- 
dependent dissipation £(i) (the latter was absent in (T^ |) 
are given by 



7(i) = ao + o>i sin(wi) 
C(t) = b[l - cos(wt)] . 



(10) 



Let e = 27rAv, / = § /„/(*)#> and /(r) = e[/(r) - 
fa], where r is defined below Eq. ([I]), for any function /. 
Applying this operation to the functions 7(er) and C( er ) 
brings Eq. (JOj) in the form of Eq. |T]) (upon subsequently 
dropping the tildes) and allows us to apply Eq. (|7|) to this 
setting. > £ cr i t rs 5.8 (in the regime of instability) 
[HI, [33] for more than half the period. Figure Q] shows 
the maximum intensity \u\ 2 of the field in the full equa- 
tion ^ in gray (green in the color online version) and 
the intensity |w| 2 of the averaged equation ([7]) in black. 
It clearly illustrates the strong correlation between the 
latter and the average (over a fast period) of the former. 
In Fig. [5J we illustrate this agreement over the initial few 
periods of the macroscopic dynamics of the equation for 
several different values of b = bo x 10 -4 , with 60 = 1, 2, 3. 
It is clear that as the magnitude of the dissipation is in- 
creased, the amplitude of the nonlinear solution sustains 
a stronger decrease and a weaker oscillatory behavior. 
Nevertheless, the agreement between the average of the 



full dynamics (over the fast scale) and the proposed aver- 
aged dynamics remains essentially satisfactory in all the 
studied cases. 



IV. CONCLUSIONS 

In conclusion, we have considered the physically real- 
istic case of periodic dissipative dynamics in the setting 
of periodically-managed nonlinearity. We have argued 
on physical grounds that the inclusion of dissipation is 
relevant for both optically layered media and for Bose- 
Einstein condensates crossing Feshbach resonances. We 
have systematically generalized earlier works by obtain- 
ing a PDE that incorporates both the average and the 
fluctuating parts of the dissipation. We showed that 
this generalized PDE model, which is valid in both one- 
dimensional and multi-dimensional settings, is in good 
agreement with the average of the original dynamical 
equation for different dissipation characteristics within 
a prototypical two-dimensional nonlinearity management 
setting proposed earlier. 

As discussed in the introduction, in addition to the 
case of temporally-dependent nonlinearity explicitly con- 
sidered here, there has been a lot of recent attention on 
spatially-dependent microstructures in the nonlinearity. 
It would thus be interesting to extend our considerations 
to the case of spatially-dependent nonlinearity prefac- 
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tors and to derive the corresponding "averaged" dynam- 
ics (provided that the nonlinearity prefactor varies over 
a fast spatial scale). Extending such averaging consid- 
erations to spatially and spatio-temporally varying non- 
linearities is an interesting endeavor under current con- 
sideration, and relevant results will be reported in future 



studies. 
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